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ABSTRACT 

Low-order, non-axisymmetric p-modes (also referred as inertial- acoustic modes) 
in hydrodynamic accretion discs around black holes are plausible candidates for high- 
frequency quasi-periodic oscillations (QPOs) observed in a number of accreting black- 
hole systems. These modes are trapped in the inner-most region of the accretion disc, 
and are subject to global instabilities due to wave absorption at the corotation reso- 
nance (where the wave pattern frequency w/m equals the disc rotation rate il), when 
the fluid vortensity, C — K^/(2r2S) (where k and S are the radial epicyclic frequency 
and disc surface density, respectively), has a positive gradient. We investigate the 
effects of disc magnetic fields on the wave absorption at corotation and the related 
wave super-reflection of the corotation barrier, and on the overstability of disc p- 
modes. In general, in the presence of magnetic fields, the p-modes have the character 
of inertial-fast magnetosonic waves in their propagation zone. For discs with a pure 
toroidal field, the corotation resonance is split into two magnetic resonances, where 
the wave frequency in the corotating frame of the fluid, uj = lj ~ niQ, matches the slow 
magnetosonic wave frequency. Significant wave energy /angular momentum absorption 
occurs at both magnetic resonances, but with opposite signs, such that one of them 
enhances the super-refiection while the other diminishes it. The combined effect of 
the two magnetic resonances is to reduce the super-reflection and the growth rate of 
the overstable p-modes. Our calculations show that even a subthermal toroidal field 
(with the magnetic pressure less than the gas pressure) may suppress the overstabil- 
ity of hydrodynamic {B = 0) p-modes. For accretion discs with mixed (toroidal and 
vertical) magnetic fields, two additional Alfven resonances appear, where w matches 
the local Alfven wave frequency. The effect of these additional resonances is to fur- 
ther reduce or diminish the growth rate of p-modes. Our results suggest that in order 
for the non-axisymmetric p-modes to be a viable candidate for the observed high- 
frequency QPOs, the disc magnetic field must be appreciably subthermal, or other 
mode excitation mechanisms are at work. 

Key words: accretion, accretion discs - hydrodynamics ~ MHD - waves ~ X-rays: 
binaries - black hole physics 

1 INTRODUCTION 

Quasi-periodic oscillations (QPOs) in black-hole (BH) X-ray binaries have been intensively studied since the launch of NASA's 
Rossi X-ray Timing Explorer (Swank 1999). The observed QPOs can be divided into two classes. The low-frequency QPOs 
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(about 0.1-50 Hz) are common, observable when the systems are in the hard state and the steep power-law state (also called 
"very high state"; see Done et al. 2007). They typically have high amplitudes and high coherence {Q <; 10), and can vary in 
frequency on short timescales (minutes). The high-frequency QPOs (HFQPOs, 40-450 Hz), however, have only been observed 
in the steep power-law state in seven BH X-ray binaries. They typically have low amplitudes (0.5 — 2% rms at 2-60 keV) 
and low coherence (Q ~ 2 — 10), with small variations in frequency (less than ~ 15%), and appear alone or in pairs with a 
frequency ratio close to 2:3 (see Remillard & McClintock 2006 for a review). QPOs in X-ray emission have also been detected 
in the active galaxy RE J1034-I-396 (Gierlinski et al. 2008; Middleton et al. 2009) as well as in several Ultra-Luminous X-ray 
sources (ULXs) (e.g., Strohmayer & Mushotzky 2009; Feng, Rao & Kaaret 2010). These are likely the "supermassive" or 
"massive" analogs of the QPOs detected in Galactic BH X-ray binaries (see Middleton & Done 2010). 

It is widely accepted that HFQPOs are associated with the dynamics of BH accretion flows near the Inner-most Stable 
Circular Orbit (ISCO). Beyond this general consensus, the physical origin of HFQPOs remains unclear, despite their well- 
defined phenomenology. In recent years, a number of models/theories with various degrees of sophistication have been proposed 
to explain HFQPOs [see Sec. 1 of Lai & Tsang 2009 (hereafter LT09) for a review of most of the existing models]. Of particular 
interest to the present paper is the relativistic diskoseismic oscillation model, in which general relativity (GR) effects produce 
trapped oscillation modes in the inner region of the disc (e.g., Okazaki et al. 1987; Nowak & Wagoner 1991; see Wagoner 
1999 and Kato 2001 for reviews). One can divide various diskoseismic modes according to their vertical structures (see Fig. 1 
of Fu & Lai 2009); The so-called p-modes (also called inertial-acoustic modes) have zero node in their wavefunctions along 
the vertical direction, while the g-modes (also known as inertial modes or inertial-gravity modes) have at least one vertical 
node^. Diskoseismic g-modes have received wide theoretical attentions because they can be trapped by the GR effect without 
relying on special inner disc boundary conditions (Okazaki et al. 1987), and they may be resonantly excited by global disc 
deformations (warp and eccentricity) through nonlinear effects (Kato 2003b, 2008; Ferreira & Ogilvie 2008). Kato (2003a) and 
Li, Goodman & Narayan (2003) (see also Zhang & Lai 2006 and Latter & Balbus 2009) showed that the non-axisymmetric 
g-mode that contains a corotation resonance (CR, where the wave patten frequency uj/m equals the background rotation 
rate O; here uj is the mode frequency and m is the azimuthal mode number) in the wave zone is heavily damped (see Tsang 
& Lai 2009a for a similar study for the c-modes). Thus the only non-axisymmetric g-modes of interest are those trapped 
around the maximum of the {Q, + K/m) profile (where k is the radial epicyclic frequency). However, the frequencies of such 
modes, cj ~ mf2(risco), are too high (by a factor of 2-3) compared to the observed values, given the measured mass and the 
estimated spin parameter of the BH (Silbergleit & Wagoner 2008; Tassev & Bertschinger 2007). Another potential problem 
with the g-modes is that their trapping zones (which arise from the maximum of the k profile near the ISCO) can be destroyed 
by a weak (sub-thermal) magnetic field (Fu & Lai 2009). For example, the cavity of the m = Q g-mode for a non-spinning 
BH vanishes when i3^/(47rpc^) <; 0.1 (where is the vertical magnetic field, p is the density and Ca is the sound speed). 
This elimination of the g-mode cavity is probably the reason for the disappearance of g-modes in recent MHD simulations 
(O'Neill, Reynolds & Miller 2009). Moreover, several numerical simulations suggest that g-modes may not survive vigorous 
disc turbulence driven by magneto-rotational instabilities (MRI) (e.g. Arras, Blaes & Turner 2006; Reynolds & Miller 2009). 

The focus of this paper is on p-modes. In an unmagnetized disc, non-axisymmetric p-modes are trapped between the 
inner disc edge (at the ISCO) and the inner Lindblad resonance (ILR) radius (where Co = uo — mQ — —k). So the existence of 
global p-modes does require a partially reflective inner disc boundary (see LT09). LT09 showed that such non-axisymmetric 
p-modes can be overstable for two reasons: (1) Because the spiral density waves inside the ILR carry negative energies while 
those outside the outer Lindblad resonance (OLR) (where iZ) = k) carry positive energies, the leakage of the p-mode energy 
when the wave tunnels through the corotation barrier (between ILR and OLR) leads to mode growth. (2) More importantly, 
p-modes can become overstable due to wave absorption at the corotation resonance (CR, where uj — 0). Near the BH, the 
radial epicyclic frequency k reaches a maximum before decreasing to zero at the ISCO. This causes a non-monotonic behavior 
in the fluid vortensity, — k^/(2ES1) (where E is the surface density of the disc), such that dC,/dr > inside the radius where 
C, peaks. It can be shown that the sign of the corotational wave absorption depends on the sign of d(/dr (Tsang & Lai 2008; 
LT09; see also Goldreich & Tremaine 1979) Thus, p-modes with positive vortensity gradient at the corotation radius can 
grow in amplitudes due to corotational wave absorption. The resulting overstable non-axisymmetric p-modes were suggested 
to offer a possible explanation for HFQPOs (LT09; Tsang & Lai 2009b). 

Realistic BH accretion disks are expected to contain appreciable magnetic fields as a result of the nonlinear development 



^ Note that other than similar mathematical forms of dispersion relations, disk g-modes have no relation to stellar g-modes which are 
driven by buoyancy. 

^ The so-called c-modes (corrugation modes) also have at least one vertical node, but they have low frequencies and are not relevant to 
HFQPOs. 

^ This applies to barotropic discs, for which the pressure is a function of the density p only. For a disc with radial stratification, the 
radial entropy gradient affects wave absorption at the CR (see Lovelace et al. 1999; Baruteau & Masset 2008), and one can define an 
effective vortensity, (^eft = C/S'^^^ (where S = P/E'" and F is the 2-dimensional adiabatic index), that plays a similar role as C (see Tsang 
& Lai 2009b). 
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of the MRI (e.g., Balbus & Hawley 1998). Starting from Hawley et al. (1995), Brandenburg et al. (1995) and Stone et al. (1996) 
in the 1990s, the property of the MRI turbulence have been studied using numerical simulations under different conditions 
and setups (local shearing box with/without vertical stratification, with/without explicit viscosity and resistivity, different 
equations of state, global simulations; e.g.. Miller & Stone 2000; Hirose et al. 2009; Guan et al. 2009; Simon et al. 2009; Davis 
et al. 2010; Promang 2010; Longaretti & Lesur 2010; Sorathia et al. 2010). Although the precise mechanism of the nonlinear 
saturation of the MRI remains unclear (sec Pcssah & Goodman 2009), the simulations generally show that turbulent discs 
are dominated by large-scale toroidal fields approaching (but less than) equipartition, with somewhat smaller poloidal fields. 
Therefore, it is important to understand how magnetic fields change the characters of the OR and the overstability of BH disc 
p-modes. These are the main issues we will address in this paper. 

We have recently shown (Fu & Lai 2009), based on local analysis of diskoseismic waves, that the propagation diagram of 
p-modes, (particularly the wave zone between inner disc edge (at nsco) and the ILR are not significantly modified by disc 
magnetic fields. Thus, we expect that the p-mode frequencies in a magnetic disc are not much different from those in a zero- 
field disc. However, the effect of magnetic field on the corotational wave absorption is more complicated. As we will see, for a 
pure toroidal field, the OR is split into two magnetic slow resonances, where the Doppler-shifted mode frequency u matches 
the frequency of slow MHD waves propagating along the field line (see Sakurai et al. 1991; Goosscns et al. 1992; Terqucm 
2003). With a mixed (poloidal and toroidal) magnetic field, two additional Alfven resonances appear, where uj matches the 
Alfven frequency associated with the toroidal field (Keppens et al. 2002). Thus, in general, the original OR is split into four 
resonances (the inner/outer magnetic slow resonances and the inner/outer Alfven resonances), which exhibit a rather rich 
and complex dynamical behavior. 

We note that Tagger and collaborators (Tagger & Pellat 1999; Varniere & Tagger 2002; Tagger & Varniere 2006) developed 
the theory of accretion-ejection instability for discs containing large-scale poloidal magnetic fields with strengths of order 
equipartition. Most of their analysis were restricted to vertical fields threading an infinitely thin disc with vacuum outside. 
They showed that the long-range behavior of the perturbed vacuum field outside the disc greatly enhances the overstability 
of non-axisymmetric disc oscillation modes. It is currently unclear whether the magnetic field configurations adopted by 
Tagger et al. can be realized by MRI turbulence - such field configuration may be produced by magnetic field advection from 
large radii (e.g. Lovelace et al. 2009). It is also unclear if the accretion-ejection instability is robust when more general disc 
field configurations are assumed. By contrast, in this paper we focus on the effects of magnetic fields inside the disc. Since 
the p-modes we are interested in have no vertical structure (i.e., kz =0), our disc model is essentially a cylinder, with the 
unperturbed disc properties (e.g., the density p, sound speed Cs, rotation Q and magnetic fields and Bz) depending only 
on r (the cylindrical radius). While the magnetic field configurations considered in this paper are prone to MRI, the p-modes 
are not directly affected by the MRI because kz = 0. 

Our paper is organized as follows. In Sec. 2, we introduce the equilibrium state of the disc and general MHD equations with 
comments on the equation structure. Wc present in Sec. 3 the effective potential of wave propagations for different magnetic 
field configurations. Section 4 examines how super-reflection across the corotation barrier is modified by the magnetic field. 
We present numerical calculations of global p-modes of magnetized BH accretion discs in Sec. 5 and conclude in Sec. 6. 



2 SETUP AND BASIC EQUATIONS 

We consider a non-self-gravitating accretion disc, satisfying the usual ideal MHD equations 

^+V.ipv) = 0, (1) 
1^ + . V)t; = -i VH -V^+J-iB- V)B, (2) 

^ = V X (« X B), (3) 
where p, P, v are the fluid density, pressure and velocity, $ is the gravitational potential, and 

H.P+|! (4) 

is the total pressure. The magnetic fleld B also satisfles the equation V • B = 0. We assume the flow obeys the barotropic 

equation of state P = P{p)- 

We adopt the cylindrical coordinates (r, (p, z) which are centered on the central BH and have the z-aods in the direction 
perpendicular to the disc plane. The unperturbed background disc has a velocity fleld v = rQ{r)^, and the magnetic fleld 
consists of both toroidal and vertical components B = B^{r)4> + Bz{r)z. The gravitational acceleration in radial direction is 
deflned as 
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so that — V"I> = —gr and g = rVt^ > 0, where SIk is the angular frequency for a test mass (the Keplerian frequency Thus 
the radial force balance equation reads 

2 dn Bl 

pg^prn^--— (6) 

To probe the dynamical properties of the magnetized flow, we apply linear perturbations to the ideal MHD Eqs. (O-® 
by assuming the perturbation of any physical variable / to have the form Sf oc e™'?'-!'^* with m being the azimuthal mode 
number and u the wave frequency. The background flow and magnetic field have no z-dependance and we will assume that 
the perturbation also has no 2;-dependance (fcz = 0). The resulting linearized perturbation equations clearly contain variables 
5v, 5p, 8P, SH and 5B. To simply the algebra, we define a new variable 

m SP B SB 

5h= — = — + — . (7) 

p p A-Kp 

Using Ai; = 5i; + ^ • Vi) = d^/dt — —icj^ + {v ■ V)^, we find that the Eulerian perturbation 5v is related to the Lagrangian 
displacement vector ^ hy Sv = —iCj^ — ri}'^r4> (prime denotes radial derivative). Also, we have (for barotropic flow) 

Sp = 5P/cl (8) 



with Cs = ^JdPjdp being the sound speed. Therefore, we can express all the perturbation quantities in terms of and 5ft, 
then further combine the perturbation equations into two equations for these two variables: 

^ = Axxir + Ax2&K (9) 

dr 

^ = AaiCr + A22&K (10) 
dr 

where 

^1 ^ ('„2 I ..2\/,-.2 ™2, .2 \lr.1 . .2\ ' I „2 i ,.2\/,~,2 , .2^ „l r.2 ™2. .2 \ ' V -I 



~ 4 2 

A ^ m 2% 
(ci +«i)(il'2 - m2a;^^)(J;2 - aj|) r'^ifij'^ - m?u\^' 

. ~2 2 2 4,{mu]\^ + CjQ.f d 2 , r 2 r>'^^rdp gdp 

A21 =iu -m UJA4, 5-^ 1- r — [ujA4, - il ) + (0JA4, - " )--r- + --r 

0J-' — m-^o;^^ dr p dr p dr 

+ 71— — 2 \/-2 ^2 2 \{~2 2T {''' [i^^ - + t^A^t^^] + - ■m'^'^M)}^ , (13) 



^ _ rCj^ [i^M ~ ^^)^^ + ^A0^^] ^ 2m{muj\^ + CbQ.) 1 dp 

'^^ ^ (ci + w^)(ii2-m2a;2^)((i2-^a;|) (ci + v^)(il)2 - aj|) r(cD2 - m^cj^ ^) p dr ' 

In the above expressions, 

H) = Lo — mO, (15) 
is the wave frequency in the co-rotating frame, 

- - J? 

is the Alfven velocity, 

VA4, 

r ry^Wp 
is the toroidal Alfven frequency, and 

=ma;A0 = fc^ (18) 



\/ci + \/ci + tJ^ ' 

is the slow magnetosonic wave frequency for k = {m/r)cf>. 

Equations ©-((TU)) describe the linear perturbations (without vertical wavenumber) in a compressible 2D flow threaded 
by a mixture of toroidal and vertical magnetic fields. In the hydrodynamic disc limit (_B — ^ 0), they reduce to the related 
equations in Goldreich & Tremaine (1979) with no external forcing potential [see also Tsang & Lai 2008 (TL08 hereafter) 



* Here "Keplerian" does not necessarily mean the gravitational potential is Newtonian. See the end of Sec. 2. 
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and LT09]. They also agree with system (14) in Blokland et aL (2005) in the fc^ — > and barotropic hmit, although their 
equations are cast on slightly different variables {611 instead of 6h). The Hameiri-Bondeson-Iacono-Bhattacharjee (HBIB) type 
of equations presented in Blokland et al. (2005) have been derived before by Hameiri (1981) and Bondeson et al. (1987) [see 
also Mikhailovskii et al. (2009); Goossens et al. (1992) (hereafter GHS92) and reference therein] for a magnetized rotating flow 
in the absence of external gravity. The equations including gravity were obtained by Keppens et al. (2002) (hereafter KCP02) 
using the same Frieman-Rotenberg technique (Frieman & Rotenberg 1960). System (5) in KCP02 is by far the most general 
formulation (with finite B^, B^, m, kz, v^, «z and adiabatic equation of state) that in principle governs all the MHD waves 
and instabilities for a compressible MHD fluid with an equilibrium flow in cylindrical geometry. In this paper, we focus on 
the effects of magnetic fields on the corotational instability of the inertial-acoustic modes (p-modes) in accretion discs, thus 
Eqs. (P))- H10|) are our main working equations. Recently, Yu & Li (2009) considered a similar system (with finite and pure 
toroidal B field) to study the Rossby wave instability and they obtained equations with similar structure as Eqs. (P))- pop . 
Our equations in the pure toroidal B field limit {Bz — ^ 0) are the same as theirs in the fcz — >■ limit. Also, the disc system 
that we study here are similar to the one considered by Terquem (2003), where instead of solving for global modes the author 
looked into the interaction between the disc and a planet and the consequences on planet migrations. 
Evidently, the coefficients in Eqs. (|11[) -(|14[1 are singular when 



~ 2 2 2 



(19) 



and 



= Lot- (20) 

In more general situations (with kz ^0), mojA4, Eq. H19|) generalizes to {k- B) / ^/AiFf), which is the Alfven frequency, and the 
slow magnetosonic frequency uos generalizes to y^5|/(ci+tJJ)(fe-i3) / ^-inp. Following Sakurai et al. (1991) and KCP02, we shall 
call p9p the Alfven Resonance (AR) and (|20p the Magnetic slow Resonance (MR). Note that these singularities/resonances 
have also been studied in the specific context of disc-planet interactions. Terquem (2003) considered a 2D disk with a pure 
toroidal magnetic field and found that a subthermal toroidal B-field is capable of stopping the inward planetary migration 
as long as the B^{r) gradient is negative enough. This analytic finding was later confirmed by 2D MHD simulation results 
(Fromang et al. 2005). On the other hand, Muto et al. (2008), by employing shearing sheet approximation, performed 3D 
(i.e., kz 7^ 0) calculations for a disk with a pure poloidal magnetic field and showed that type I planetary migration could be 
outward if the B-filed is strong enough (superthermal) . 

In the special case of kz = Bz = 0, Eq. (|19p is no longer a real singularity of the differential equations. In this case, one 
can show, through some subtle mathematical manipulations, that the terms cj^ — m'^uj\^ in An, A12, A21 and j422 all get 
canceled. These four coefficients then reduce to 



^11 



ci + V 



A4, ' 



+ 2mil!Q 



^ ( n2 

^A<i> 



9) 



A12 = 



- 2 2 2 /2 

Lj — m c^/r 



"A^ 



+ ci 



I21 = ti; — m L0A4, 



dln{pr/Bl) v% 



dlnr 



ci + V 



A4> 



(21) 



(22) 



mvA^ 



(Sl^ - g/r + 2c? /r") 



-'Act, 



2 

+ AmCjQ.^.^^^^i^'' - g/r + 2c? /r') ), 



ci -\- V 



Atf, 



(23) 



r ci+ 



2 

2CjQ + ,"^"^1 {Q,'^ - g/r + 2c?/r^) 



ci + V 



A4, 



where the radial epicyclic frequency k is given by 



2Q d 2 
r dr 



1/2 



d\np 
rflnr 



+ 



c| + V 



-(rfi - g ~2vA4,/r), 



Acj, 



(24) 



(25) 



This peculiar disappearance of Alfven singularity can also be verified from a different perspective. GHS92 derived the jump 
conditions for and SU (called P' in their paper) across the Alfven resonance 

[Cr] oc gB = ruBz/r ~ kzB^, 

[SU] (X Bz 



There is a typo in their Eq. (20), where in the last line, the term (BgU) + Fvg) should be (B| 
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(see their Eqs. [42]- [43]). We see that when = = 0, both and SIl are continuous across the Alfven resonance, thus 
the Alfven singularity disappears 0. In a non-magnetic disc, Eqs. (|19|l - (|20p reduce to the corotation singularity (TL08). 

In Sec. 3 and Sec. 4, we will simply employ Newtonian gravitational potential, which has free-particle orbital frequency 
flK oc r~^^^ and epicyclic frequency k = Qk- This way a dimension-free analysis can be easily made without missing the 
essential physics. However, in Sec. 5, in order to make a direct comparison with results in LT09, we will use the Paczynski-Wiita 
pseudo-Newtonian potential (Paczynski & Wiita 1980) 

<f = -^, (26) 

r ~rs 

where M is the central BH mass and = 2GM/(? is the Schwarzschild radius. The corresponding free-particle orbital 
frequency and radial epicyclic frequency are 

(Id^y'^ IGM 1 
\rarl V r r — Ts 



K = nKJ^-^. (28) 

V r - rs 

Note that peaks at r = (2 -I- \/3)''s and decreases to zero at r = 3rs. Throughout this paper, we will consider thin discs 
with the sound speed Cs <^ tQk and magnetic field satisfying _B^/(47rp) ^ (rf^K)^, so that Q ~ Qk- 



3 EFFECTIVE POTENTIAL AND WAVE PROPAGATION DIAGRAM 

We can combine Eqs. (f9)l- (|10p into a single second-order differential equation on 5h 

—Sh + Ci{r)—5h + Co{r)5h = 0, 
dr^ dr 

where 
Ci{r) 



-An 



1 dAai 
A21 dr 



CQ{r) = A11A22 - A12A21 ; h 



dr 



A21 dr 



We introduce a new variable 
1 



r/ = exp 



rCi(r)dr 



Sh, 



so that Eq. (|29p can be rewritten as a wave equation for 77: 
where 

V(.) = iCi(r)^ + i^^-Co(r) 

is the effective potential for wave propagation. Considering a local plane wave solution 



rj oc exp 



kr{s)di 



Eq. (I33j then yields 
fc^ = -V{r). 



(29) 

(30) 
(31) 

(32) 

(33) 
(34) 
(35) 
(36) 



Formally, wave can propagate only in the region with V{r) < 0. Our derivation above is quite general and can be applied to 
complicated systems as the details have all been captured in the four coefficients Aij. Figures [T][3] show the effective potential 
for various situations discussed below. 



^ Although the analysis in GHS92 does not consider external gravity, one can show that the effect of gravity can be easily included 
through replacing the term pv^/r in their Eq. (5) by pv^/r + pg. This modification only changes the relations between the background 
variables {B^, v^, Cs, etc.) The jump conditions for and P' would remain the same. 



Corotational instability in magnetized discs 7 



3.1 Hydrodynamical discs 

For a non-magnetic disc, the four coefficients Aij simpfify to 
rfln(rp) 2mn 

All = 1 —, (•31) 

ar ruj 

Ai2^-^ + ^, (38) 
A21 = -D, (39) 

A22 = — —, (40 

ruj 

2 



where D — k, — Hi . Thus Co(r) and Ci(r) now read 
D m'^ 2mQ, f d Qp 



Co(r) = --- — ^ In-^ , (4f) 

Ci(r) = --f In-. (42) 
dr rp 

The effective potential then takes the following form 

Vir) = ^ + 4+'-^(^ln^)+ S^/^ (43) 
cj ru! \dr D J dr-' 

where S — D/(rp). The effective potential V{r) is singular at both the corotation resonance (CR), where a) = 0, and the 
Lindblad resonances (LRs), where D = Q (see Figs. lH3"|l . However, the singularities at the LRs are spurious singularities of the 
wave equation (i.e. they can be removed by rewriting the equation in alternative variables [e.g., Narayan et al. f987; TL08]) 
and only the CR is a non-removable singularity. The behavior of V{r) around the CR depends on dC^/dr (evaluated at the 
corotation radius rc), the gradient of the background fluid vortensity, defined by 

^-4rp 

For d(/dr < a narrow "Rossby wave zone" lies just inside Tc, while for d(/dr > the Rossby zone lies just outside rc; in 
the special case of d(/dr = 0, the corotation singularity disappears (see TL08). 



3.2 Discs with pure toroidal magnetic fields 

As mentioned in Sec. 2, for a disc with a pure toroidal magnetic field, the CR splits into two magnetic slow resonances (MRs). 
This can be seen clearly from the wave propagation diagram shown in Figs. [l][3l For B = 0, the effective potentials reduce to 
those shown in TL08. In the upper panels of Figs. [T][3l we see that for small magnetic fields (with lj Att> ^ fi), the effective 
potential profile is not affected much by the B-field except near the corotation. The blow-ups in the bottom panels show that 
the toroidal magnetic field splits the CR into two MRs, the inner magnetic resonance (IMR) and outer magnetic resonance 
(OMR), one on each side of the CR. A larger toroidal field (see the long-dashed lines) gives a wider radial separation between 
the MRs and a shallower potential well around rc. Note that independent of the sign of (d("/dr)c, there always exists an 
effective wave zone centered around corotation point and bounded by the two MRs. 

Far away from the resonances (the regions on the left and right sides of the upper panels of Figs.[T][31), the WKB dispersion 
relation 1(351) reduces to (see Terquem 2003; Fu & Lai 2009) 

^b"" ^ K,"" + kf.{ci + v%). (45) 

This describes fast magnetosonic waves modified by the disc rotation. Thus the spiral density waves (inertial-acoustic waves) 
inside the ILR and outside the OLR are inertial-fast magnetosonic waves. 

The behavior of the effective potential around the MRs is more difficult to analyze. The IMR and OMR are separated 
from the CR by the distance 

'-'-"-^^^i^)'^' ^^^^ 

(the upper sign is for the OMR and the lower sign for the IMR). Because of the complexity of the equations involved, no simple 
analytical expression for V{r) can be derived, even in the asymptotic limit (e.g., to — >■ ztujs and VAt/>/rQ ^ 1). Nevertheless, 
a careful examination of Eqs. (|2f p - (|24| l. (|30|) -(|3I | I. (|34|) and numerical calculation of V(r) show that in the neighborhood of 
the MRs, the effective potential can be written in the following form: 
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Figure 1. Wave propagation diagram in accretion discs when the background fluid vortensity gradient (evaluated at the corotation 
resonance radius) {dQ/dr)c = 0. In the upper panel, the solid line shows the effective potential V{r) as a function of r for B = discs, 
and the dashed line shows V{r) for discs with finite toroidal magnetic fields. Waves can propagate only in the region where V{r) < 0. 
The labels ILR, OLR and CR denote the inner Lindblad resonance, outer Lindblad resonance and corotation resonance, respectively. The 
lower panel shows the blow-up of the region near the CR (located at rc = 1), and the long-dashed line shows the case with a stronger 
than the short-dashed line. The labels IMR and OMR denote the inner magnetic (slow) resonance and outer magnetic resonance, 
respectively. Note that the vertical scales in the upper and bottom panels differ by a large factor. A Newtonian potential is used in 
calculating V{r). (A color version of this figure is available in the online journal.) 



V{r) = fo + ^J^ + -J^ + -J^, (47) 

where /o, /i, /2± vary slowly around the CR/MRs region. Thus, the effective potential is dominated by second-order singu- 
larities at the MRs, i.e., V{r) cc (r — rmr)~^ (see Figs. [l][3]). The coefficients /2± can be calculated numerically, and some 
examples of how f2+ and /2_ depend the magnetic field strength are shown in Fig. ID We see that in dimensionless units 
(with G — M = Tc — 1, m = 2 and Cs <C rQ, va</> ^ rQ,), both f2+ and /2_ are negative and approximately equal to —2. 
Thus, near one of the MRs (e.g., the OMR), we have fcj: = -V{r) ~ 2/{Lj-iUsf ~ 0.2/{r - Vmi)^ , or kr ~ OA/\r-rmi\. Since 
is smaller than \dkr/dr\, the WKB analysis is not really valid near the MRs. On the other hand, close to the CR radius 
(r = rc), we have V{r) ~ —i/cu^, which implies kr ~ 2uj~^. Since the separation between the IMR and OMR is 2ujs/-i^ (for 
m — 2), we find that the WKB phase variation across the CR region (but avoiding the MRs) is of order unity. Therefore, 
despite the deep effective potential at the MRs, physically there is not much a wave zone around the MRs/CR region. This 
feature is borne out in our numerical calculations of global wave modes presented in Sec. 5. 

It is worth noting with the second-order singularities associated with the MRs, the _B — >■ limit is approached in a 
non-trivial way. From Fig. [4l we see that as the magnetic field decreases, /2+ and /2- reach the same values (~ —2.25). 
Thus the singularities associated with f2+/{0j — cjs)^ and /2_/(iI; + ujs)'^ simply become as B — >■ 0. On the other hand, 

/i can be written as /i = fio'jj'^ + /iiii + /i2<i^ + ■ ■ with fu being near constant around the MRs. As B — >■ 0, we have 
/if) + /2+ + /2- = 0, and the second-order singularity term (oc l/ti^) in V{r) disappears. What is left is the term /ii/cD, a 
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first-order singularity associated witii the CR in a non-magnetic disc. Evidently, as B — )■ 0, /n is proportional to the vortensity 
gradient, d(/dr. 



3.3 Discs with pure vertical magnetic fields 

Before examining the case of mixed magnetic fields, it is useful to consider discs with pure vertical fields. Here we assume 
that Bz is constant (independent of r) to simplify the analysis. The coefficients Aij of the wave equations (f9ll- (|10p are 

cj dlnp 2m!:; 1 

^ ~ „2 I „.2 — Z ZFT^Z' (48) 



"A 



r 



^-=^^-^-7^^f^V. (50) 



"A 



v\ dlnp 2mQ, 
ci + dr ru) 

We then have 

^ D 2mD. f d , Q.p^ 
Co = — — 5 5 — -r in 



dr A 
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Figure 3. Same as Fig.[T] except for the case {dC,/dr)c > 0. 



clv\ / d\ap \ v\ / dlnp \ / dln(yl2i/r) \ ^ d f v\ dlnp \ 



(,2 _|_ ^2 ^ y j,2 ^ ^2 ^ y y y dr \ + dr 

where 

A = (c? - 4)/(c^ + 4) (54) 

and D = — Cj^ . The resulting effective potential reads 

= ^ ^ ^ / d ^7/ X 5V2^5-V2 

ci~\-v\ ro; \dr A21 / dr^ 

cWa / dlnp y ^ ^jj / dlnp X / dln(^2i/r) \ d / t;^ dlnp X 
r^il)^ c1+v\\dr) ci+v\\drj\ dr J dr + v'\ dr J ' 

where S — A21/ (rp). Compared with the V{r) in the 5 = case, we see that a constant vertical field does not split the CR. 
However, in contrast to the first-order corotation singularity in the B — case, the CR in discs with _Bz 7^ manifests 

mainly as a second-order singularity (~ l/cj^). 



3.4 Discs with mixed magnetic fields 

For a disc with mixed magnetic fields, the expression for the effective potential V(r) is more complicated than the pure toroidal 
field case. In Fig. O we show a sketch of V{r) when both Bz and B^ are non-zero. We see that in addition to the two MRs 
(IMR and OMR) (already present in the pure toroidal field case), two Alfven resonances (ARs) appear when d) — ±muiA^- 
Since mujA(/> > ^a, these ARs lie farther away from the corotation radius than the MRs, at the distance 
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Figure 4. Numerical values of /2+ and /2— [see Eq. I I47II ] as a function of the dimensionless ratio D40/(rf2) (evaluated at r = Tc), 
which specifies the toroidal magnetic field strength in the disc. Dimensionless units are adopted, where G = M = rc = 1, m = 2, and 
Cs/(rO) <C 1. Different line types represent different disc density and magnetic field profiles. 



Unlike the MRs, the ARs are first-order singularities of the wave equation, i.e., V{r) oc — mujA(f,) oc l/(r — rar) near the 
ARs. 



4 SUPER-REFLECTION OF THE COROTATION BARRIER 

As shown in Sec. 3, away from the CR/MRs/ARs region, the WKB dispersion relation is simply 

Cj'^ = K,^ + ki{(?s + va). (57) 

Thus waves can propagate either inside the ILR (where Cj — — k) or outside the OLR (where Cj — k), and the region between 
the ILR and OLR represents the "corotation barrier", where waves are evanescent except for various possible singularities 
associated with the CR/MRs/ARs. We will be interested in the wave modes trapped between the disc inner edge and ILR. As 
in the hydrodynamic {B — 0) case, the growth rates of such p- modes depend directly on the wave reflectivity of the corotation 
barrier (TL08, LT09). 

Consider launching a density wave Sh oc exp(— i krdr) from a radius inside the ILR (r < tilr) (assuming kr > 0; 
the minus sign arises from the fact that group velocity and phase velocity have different signs where r < riLR, see TL08). 
The wave carries negative energy since its pattern speed oj/m is smaller than background rotation rate Q,. When the wave 
impinges on the corotation barrier, it is partly bounced back as a reflection wave Shn oc 7?.exp(i J ^ krdr), which carries 
negative energy, and partly transmitted through the evanescent zone and travels beyond the OLR as a transmission wave 
(5/iT oc Texp(i krdr), which carries positive energy. TZ and T are reflection and transmission coefficients, respectively (see 
TL08). Let the net absorption of wave energy at the CR, MRs and ARs be Oabs- Energy conservation then implies 

17^f = i + |r|' + OabB. (58) 
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Figure 5. Effective potential around the corotation/magnetic resonance region for a disc with mixed magnetic fields (solid line). The 
dotted line shows V{r) for the case with a pure toroidal field. The corotation resonance (CR), inner/outer magnetic (slow) resonances 
(IMR and OMR) and inner/outer Alfvcn resonances (lAR and OAR) are indicated. Note that the MRs represent second-order singularities 
in V(r), while the ARs (which exist only for the mixed field case) are first-order singularities. 



Obviously, wave transmission through the corotation barrier (the |T|^ term) always tends to increase |7?.|^ beyond unity. 
However, the resonant wave absorption term T>ahs can be either positive or negative. 



4.1 Non-magnetic Discs 

TL08 showed that, for a. B — Q barotropic disc, the wave super-reflection (|7?.p 
absorption, which depends on the dimensionless parameter 



pK ar J 6 



1) is dominated by the corotational wave 



(59) 



where is the vortensity [see Eq. (|44|l ]. /3, p and a are defined by /3 = Cs/[rQ.), Q, oc and p oc r~'^ (around the CR), 
respectively, and the subscript "c" implies that the quantity is evaluated at r = r^. The second equality in Eq. (|59|) assumes 
(Newtonian) Keplerian discs (p = 3/2). Importantly, the corotational wave absorption is proportional to u or d(/dr. This can 
be understood from the fact that for v > {v < 0), the Rossby zone lies outside (inside) Vc (see Figs. IH3")) . implying a positive 
(negative) wave energy absorption. Thus, if neglecting |Tp, we have 



1711 



1 oc oc 



(60) 



Tsang & Lai (2009b) further generalizes the analysis to non-barotropic flows with radial stratification, in which case an 
effective vortensity (which depends on the radial entropy profile of the disc) plays a similar role as (. 
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Figure 6. Reflection coefHcient as a function of u for a Newtonian Keplerian disk. Different lines represent different magnetic field 
strengths, measured by the dimensionless parameters f)^ = {vj^^/rQ)c and bz = {vjHz/rCl)c. The other parameters are m = 2, oc 
Bz oc r"? with 5 = and /3 = c,,/(rf!) = 0.1. The parameter u is defined according to Eq. II59I I. 



4.2 Numerical calculations of reflectivity 



For discs with magnetic fields, we are not able to derive analytical expressions for the reflectivity. Thus we calculate TZ 
numerically. 

For definiteness, we consider a (Newtonian) Keplerian disc with 



O, = K oc r 



-3/2 



peer Ba, cc Bz cc r'', = /3, 



(61) 



where a, q and /3 are constants. The strength of the field is characterized by the ratio — VA4>/{'rQ.) and = VAz/irQ.), 
evaluated at rc. We choose the wave frequency uj = uir + ii^i, with Q < oji <Si o^r- The transmitted wave outside the OLR takes 
the form 



5h oc A exp l^i j fc,. drj , 

where fcr > is given by Eq. (|36() or approximately by Eq. (|45p . and A{r) can be obtained from Eq. \'i2\ with i] oc k, 
Thus the boundary condition at rout > '"olr is 

dSh i 
dr \r. 



: Sh {ikr + A' /A) 



(62) 



-1/2 



(63) 



where A' = dA/dr. Starting from r = rout, we integrate Eqs. (f9|l- (|10p inwards to a radius inside the ILR. At r = rin < hlr, 
the wavefunction has the form 



cxp ( — i / krdr] +71 exp ( i / krdr 



Shoe A 

The reflectivity TZ can then be extracted from the equation 

_ (-ifc + A'/A)Sh - Sh' 
' ' ^ {ik + A'/A)Sh-Sh' ^^^^ ■ 

In practice, at rin (rout) sufficiently inside riLR (outside roLR,), the A' /A term can be dropped since kr ^ j^'/^j. 



(64) 



(65) 
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Figure [6] shows the numerical results of the reflection coefficient as a function of u, defined by Eq. (|59|) . For B = 0, 
our result recovers that obtained by TL08. We see even a relatively weak (sub-thermal, with VA4, ^ Cs) toroidal field can 
significantly affect the reflectivity. As the magnetic field increases, the reflectivity decreases regardless of the sign of u. A 
finite vertical field has a similar effect, although it is sub-dominant compared to the toroidal field effect. We will show in the 
following section that this result is consistent with the stability property of global disc p-modes. 



5 GLOBAL P-MODES OF BLACK-HOLE ACCRETION DISCS 

In this section, we present our calculations of the global non-axisymmetric p-modes in BH accretion discs. Recall that for a 
nonmagnetic disc, these are the inertial-acoustic modes partially trapped between the inner disc edge (at nn = risco) and 
the ILR. For a disc with finite magnetic fields, these become inertial-fast magnetosonic waves. 

We employ the standard shooting method (Press et al. 1992) to solve Eqs. (f9|- (|10() for the complex eigenfrequency 
uj = ujr + The boundary condition at rout > ^olr is the outgoing-radiative condition discussed in Sec. 14.21 [see Eq. (|63p ]. 
At the inner boundary rin = risco, we apply a free surface boundary condition, 

= An/p = 0, (66) 

i.e., the Lagrangian perturbation of the total pressure vanishes. This boundary condition corresponds to a idealized situation 
in which no wave energy is lost at the inner disc radius. As discussed in LT09 (and references therein), the inner edge of a 
real BH disc is more complicated and may involve wave energy loss due to radial infall of the accreting gas. Since our focus 
here is to understand how disc magnetic fields affect the corotational instability, we will adopt the simple boundary condition 
H66|) for all the numerical calculations presented in this section. 

As mentioned in Sec. 2, we use the Paczynski-Wiita pseudo-Newtonian potential [Eq. (|26|) ] to mimic the GR effect. The 
disk density profile and magnetic field profile take the power-law forms 

pocr"", B^ocBzOcr''. (67) 
The thermal effect and the magnetic field effect are characterized by the dimensionless parameters 

with /3 being constant. 

The angular momentum flux carried by the wave across the magnetized disc is given by (e.g. Pessah et al. 2006) 

F{r) = Tvr'^pRe{Svr5v;) - K'^Re{SBr5B;), (69) 



where * denotes complex conjugate. The first and second terms in (|69[) are related to the Reynolds stress and Maxwell stress, 

respectively. The solution of Eqs. (f9))- (|10p gives (,r{r) and Sh{r). The expressions for Svr, Sv^, SBr and 5B^ can be derived 
from the original linearized perturbation equations (HJ-Q: 

5vr = —iu^r, (70) 

S,, = + ^A,^ ^r+(^ + ^Ar.) Sh, (71) 
\ 2il LJ ruj I \rLj ru! I 

SB. = (72) 
r 

5B4, = i-q^ - ^iiB*.) Cr - Ai2B^Sh. (73) 



5.1 Results for discs with pure toroidal magnetic fields 

Figure [57H gives an example of the eigenfunctions for overstable p-modes in magnetized BH accretion discs. For comparison, 
the corresponding results for a i? = disc are also shown. Clearly, with a small B field (6^ <^ 1), the wavefunctions away 
from the corotation/magnetic resonances (CR/MRs) are only slightly modified. The wave action is still trapped in the inner 
disc region. The most noticeable effects occur near the MRs, where both 5h and Svr experience significant variations. In 
this example, the 5 = disc has dC,/dr > at the CR, so that corotational wave absorption leads to a mode growth rate 
uji ~ 0.0027nin. The inclusion of a small toroidal B field, fe^ = 0.01 [corresponding to {vA<t,/r^)t^ = 0.0175], however, reduces 
the mode growth rate to Ui ~ O.OOlSr^in. Note that the disc sound speed in the example is Cs — O.lrQ, so the 6^ = 0.01 
corresponds to {vAcj>/cs)c = 0.175. 

The origin of the diminished mode growth rate due to magnetic fields can be understood by examining the angular 
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Figure 7. Example wavefunctions of m = 2 p-modes in BH accretion discs. The disc density profile is p oc (i-e., cr = 1) and the sound 
speed is Cs = O.lrQ (i.e., /3 = 0.1). The left columns are for a, B = disc, with the eigenvalue uj = (0.9324 + i0.0027)nin (where 0;^ is the 
disc rotation rate at r = rj,! = ngco); the right panels are for a disc with = (fyi</)/''^)in = 0-01 [corresponding to {vacI>/^^)c = 0.0175] 
(with Bif, independent of r, i.e., q = 0), with u} = (0.9312 + j0.0018)nin. In each subfigure, the upper panels show different variables as a 
function of r in the full fluid zone. The solid and dotted lines depict the real and imaginary parts, respectively. The bottom panels zoom 
in the region near the corotation resonance, with the vertical lines indicating the locations of the CR (left) and magnetic slow resonances 
(right). The vertical scales of the wavefunctions are arbitrary, with 5h{rout) = 1. 



momentum flux carried by the wave modes. Figure [8] gives some examples. We see that in the B = case, there is a jump of 
the angular momentum flux across the CR, 

AF = F(rc-) - F(rc+), {B = 0). (74) 

The fact that AF > implies a positive wave energy (angular momentum) absorbed at the CR [i.e., ©abs > in Eq. (|58|) ]: 
this leads to super-reflection and is the main driver for the overstability of hydrodynamic p-modes. For discs with a finite B^, 
the CR is split into two MRs (the IMR and OMR), and we see from Fig. [8] that significant flux jumps occur at both the IMR 
and OMR. The flux jump at the IMR, AFimr ~ F(nMR-) — i^(nMR+), is negative, implying a negative angular momentum 
absorption. The flux jump at the OMR, AFomr = J'(roMR-) — F(roMi{.+ ), however, is positive, implying a positive angular 
momentum absorption. The net angular momentum jump across the MRs region is then 
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Figure 8. Angular momentum flux carried by the wave as a function of r. The disc parameters are a = 1, q = and /3 = 0.1, and the wave 
modes have m = 2. Panel (a) shows the result of a iJ = disc [with the mode frequency u) = (0.9324 + i0.0027)r2in] with the vertical line 
indicating the corotation resonance. Panel (b) is for a magnetized disc with = 0.01 [corresponding to {vA^/rQ)c = 0.0175] [left column; 
with ui = (0.9312 + i0.0018)nin] or = 0.015 [corresponding to {vA4,/rO.)c = 0.026)] [right column; with ui = (0.93 + i0.00065)r2in]. 
The solid and dotted curves show the total angular momentum flux and the flux carried by the fluid motion only (the Reynolds stress), 
respectively. The bottom panels are the blow-up of panel (b) near the magnetic resonances whose locations are indicated by the two 
vertical lines. 



AF = F(riMR-) - F(roMR+) = AFimr + AFqmr. (75) 

Note that although both [AFimrI and | AFomrI are much larger (by orders of magnitude) than AF in the B = Q case, the net 
jump AF is smaller for discs with finite B^. As B^ increases, |AFimr| and [AFomr! both decrease, and AF becomes smaller 
and may even change sign (See Fig. [9|. This decrease of wave absorption at the MRs due to finite magnetic fields results in 
a reduced or even diminished super-reflection (see Sec. 4), leading to a reduction of the growth rate and even stabilization of 
disc p-modes. 

Figure [9] shows the p-mode growth rate as a function of and the corresponding net angular momentum jump across 
the CR/MRs region. Obviously, the decrease of AF with increasing directly correlates with the decrease of the mode 
growth rate. Note that when AF is zero or slightly negative, the p-mode still has a finite growth rate. This is because wave 
transmission across the corotation barrier always tends to increase the reflectivity [see Eq. (|58l) ]. thereby promoting the the 
mode growth. With a further increase of b^, the mode growth becomes completely suppressed. For example, for discs with 
Ca = O.lrfi, the angular momentum absorption at the CR/MRs change sign at b^ = 0.0132, while the mode growth rate 
becomes zero at ~ 0.017 , or («A0/cs)iii — 0.17 [corresponding to («A(^/cs)c = 0.3]. For discs with Cs = O.OSril, the mode 
growth rate vanishes at b^ ~ 0.01, or (tJA0/cs)in — 0.2 [corresponding to (va4>/cs)c = 0.34]. 

Figure fTOl shows the real and imaginary parts of the m = 2 mode frequency as a function of the dimensionless field strength 
va4>/('i'^) (evaluated at Tc) for different disc sound speeds and different values of q (measuring the slope of B^). Clearly, for 
such a small field strength, the real mode frequency is approximately independent of B^ and q, since the propagation zone of 
the wave mode is hardly modified by the magnetic field. The mode growth rate, however, is significantly affected because of 
the modification to the CR. Not surprisingly, Im(aj) depends on B^ mainly through the ratio {va4>/t^)c- 



5.2 Results for discs with mixed magnetic fields 

The effects of vertical magnetic fields and mixed fields on the disc p-mode growth rate are illustrated in Fig. [TT] and Fig. 1121 
We see that the real mode frequency depends very weakly on the field strength unless va becomes comparable to or larger 
than Cs . The solid line in Fig. Illb shows that for the pure vertical field case, the p-mode growth rate decreases with increasing 
_Bz. This is consistent with our finding in Sec. 4.2 that a finite B^ tends to reduce the reflectivity (see Fig. [6]). When the 
toroidal field is also included, we see the similar trend (see the dotted and dashed lines). Thus both vertical and toroidal 
magnetic fields tend to suppress the corotational instability. The effect of the toroidal magnetic field is somewhat larger that 
the vertical field, given that the required toroidal field strength B^ to completely suppress the mode growth is smaller than 
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Figure 9. The net angular momentum flux jump across the two magnetic resonances (upper panel) and the mode growth rate (lower 
panel) as a function of the dimensionless toroidal magnetic field strength 6^ = {vj^^/rQ)i^. The solid lines are for disc sound speed 
Cg = O.lrfi, and the dashed line in the lower panel is for Cg = O.OSrn. The other parameters are the same as in Figs . [STTI and 151 m = 2, 
a=l, q = 0. 



the required B^. This is understandable since the p- modes have no vertical structure (with =0), and fluid motion in the 
disc does not directly bend the vertical held lines. 

In Fig. 1121 we compare the angular momentum flux profiles for the pure toroidal field case and the mixed field case. From 
the upper panel of this figure we see that for the mixed field case, the net fiux jump across the resonance zone is smaller than 
the pure toroidal field case. This explains the somewhat smaller mode growth rate (oji/fiin = 0.0017) for the mixed field case 
as compared to the pure toroidal field case (cJi/Slin = 0.0018). Note that the real part of the eigenfrequency is also modified 
slightly by the vertical field. This gives rise to the shift of "plateau" in the bottom panel of Fig. 1121 



6 CONCLUSION 

In this paper we have investigated the effects of magnetic fields on the corotational instability of the p-modes (also called 
inertial-acoustic modes for the B — case) in BH accretion discs. These modes have no vertical structure, are trapped between 
the inner disc edge and the inner Lindblad resonance, and have the character of inertial-fast magnetosonic waves in their 
propagation zone. Previous works have shown that in a hydrodynamic disc (with B = 0) these modes can be overstable due to 
wave absorption at the corotation resonance (CR) when the vortensity (or its generalization for non-barotropic flows) of the 
disc has a positive gradient (TL08, LT09, Tsang & Lai 2009b). In this paper we found that the inclusion of a finite toroidal 
field B^ splits the CR into two magnetic (slow) resonances (MRs), where the wave frequency in the rotating frame of the 
fiuid, uj — mfl, matches the slow magnetosonic wave frequency. At the inner and outer MRs, the angular momentum fiux 
carried by the wave undergoes significant change, indicating angular momentum (and energy) absorption at both MRs. One 
of the resonances absorbs positive angular momentum while the other emits positive angular momentum. Independent of the 
background flow vortensity gradient, the net angular momentum absorption across the resonance region is always reduced 
or becomes more negative in discs with B,j, 7^ compared to the B — case. This leads to a reduced growth rate of the 
p-modes. Our calculations showed that the hydrodynamically overstable inertial-acoustic modes can be completely stabilized 
by the toroidal fleld at relatively small field strengths (see Fig. llOp For example, for discs with sound speed Cs = O.lrfi, the 
m — 2 p-mode growth rate vanishes at VAci>/{r^) ~ 0.03 (evaluated at the corotation radius rc; where VA4, ~ B^j^f^Kp is 
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Figure 10. Frequencies of the m = 2 p-modes as a function of the dimensionless toroidal magnetic field strength, {va<i>/i'Q,)c, with 
the upper and lower panels showing the real and imaginary parts, respectively. Different curves correspond to different values of g = 
d\n B^/ d\nr (i.e., the slope of profile) and sound speed c^. The disc density profile is p oc r~^, i.e. a = I. 



the Alfven speed associated with B^), corresponding to {vA<t,/cs)c ~ 0.3. We also carried out numerical calculations of the 
reflectivity of a wave approaching the corotation from small radii. We showed that the reflectivity is reduced by the toroidal 
magnetic field, in agreement with our global mode calculation results. 

For discs with a pure vertical field, the corotation resonance persists for the p-modes, with no additional magnetic 
resonances. This is understandable since p-modes have no vertical structure (fcz = 0). With no bending of the field lines, 
the vertical field simply modifies the background pressure, or equivalently, changes the effective background sound speed. 
However, when 7^ 0, the effective potential for wave propagation contains a second-order singularity term, in addition to 
the first-order singularity already present in a = disc. 

For discs containing mixed (toroidal and vertical) magnetic fields, the corotation resonance is split into four resonances: 
in addition to the inner/outer MRs (already present in discs with pure toroidal fields), two Alfven resonances appear, where 
Lj — mQ matches the local Alfven wave frequency. We showed that the effect of these additional resonances is to further reduce 
the super-refiection and the growth rate of the disc p-modes (see Fig. Overall, the toroidal field has a larger effect on the 
corotational wave absorption and the mode growth rate than the vertical field. 

The main finding of our paper is that overstable non-axisymmetric p-modes tend to be stabilized by disc magnetic fields. 
This implies that, in order for p-modes to explain HFQPOs observed in black-hole X-ray binaries, the disc magnetic field 
must be sufficiently weak (with the magnetic pressure appreciably less than the thermal pressure) so that the corotational 
instability is not completely suppressed. Alternatively, some other mode excitation mechanisms are needed. One possibility 
could be mode excitations by turbulent viscosity (see Kato 2001 for a general discussion of viscous driving of disc oscillation 
modes) or by vorticity perturbations associated with the turbulence (see Heinemann & Papaloizou 2009). Another possibility 
could be the accretion-ejection instability studied by Tagger et al. (see Tagger & Pellat 1999; Varniere & Tagger 2002; Tagger & 
Varniere 2006), which involves large-scale equipartition vertical magnetic fields threading a thin disc embedded in a vacuum or 
a tenuous corona. Whether these other possible excitation mechanisms can compete with the stabilizing effects of the internal 
disc magnetic fields remains to be studied. Nevertheless, our finding is limited by the fact that the disc model we considered 
is essentially an infinite cylinder without vertical structure, whereas in reality disc background (density, magnetic field, etc.) 
varies in vertical direction. More realistic disc models certainly need to be investigated in future studies. 

Finally, we note that although in this paper we have studied the role of magnetic fields in the context of black-hole 
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Figure 11. Frequencies of the m = 2 p-modes as a function of ttie dimensionless vertical magnetic field strengtli, {v^^^ /rQ)c, with the 
upper and lower panels showing the real and imaginary parts, respectively. Different curves correspond to different values of the disc 
toroidal fields. Note that the real mode frequency depends weakly on the magnetic field (for the range of field strength considered), thus 
= 0.01 corresponds to (i)40/rf2)c ~ 0.0175 and = 0.015 corresponds to (i)^0/rf2)c ~ 0.026. The other disc parameters are: cr = 1, 
q = and Cs = O.lrSl. 



diskoseismic oscillations and QPOs, the dynamical effects of magnetic fields on the corotation resonance examined here are 
quite general. Our finding in this paper suggests that dynamical instabilities in other rotating astrophysical flows where 
the corotation resonance plays an important role may be significantly affected by magnetic fields. Examples include the 
Papaloizou-Pringle instability in accretion tori and the global rotational instability in differentially rotating stars. We will 
study some of these issues in separate papers (Fu & Lai 2010). 
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